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ABSTRACT 

Aims. This work aims at investigating the molecular gas in the surroundings of the ultra-compact HII region 
G045.47+0.05 looking for evidence of molecular outflows. 

Methods. We carried out observations towards a region of 2'x2' centered at RA=19 h 14 m 25.6 s , 
dec.=+ll°09'27.6"(J2000) using the Atacama Submillimeter Telescope Experiment (ASTE; Chile) in the 12 CO 
J=3-2, 13 CO J=3-2, HCO + J=4-3 and CS 3=7-6 lines with an angular resolution of 22". We complement these 
observations with public infrared data. 

Results. We characterize the physical parameters of the molecular clump where G045. 47+0.0 is embedded. The 
detection of the CS 3=7-6 line emission in the region reveals that the ultra-compact HII region G045. 47+0.0 has not 
completely disrupted the dense gas where it was born. The HCO + abundance observed towards G045. 47+0.0 suggests 
the presence of molecular outflow activity in the region. From the analysis of the 12 CO J=3-2 transition we report the 
presence of bipolar molecular outflows with a total mass of about 300 M©. We derive a dynamical time (flow's age) 
of about 10 5 yr for the outflow gas, in agreement with the presence of an ultra-compact HII region. We identify the 
source 2MASS 19142564+1109283 as the massive protostar candidate to drive the molecular outflows. Based on the 
analysis of its spectral energy distribution we infer that it is an early B-type star of about 15 M©. The results of this 
work support the scenario where the formation of massive stars, at least up to early B-type stars, is similar to that of 
low mass stars. 
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1. Introduction 

The formation of high-mass stars remains one of the most 
significant unsolved problems in astrophysics. Despite the 
importance that massive stars have in the structure and dy- 
namic of the Galaxy, the physical processes involved in their 
formation are less understood than those of their low-mass 
counterpart. Observationally, the main problems arise from 
the fact that they are heavily obscured by dust, are rare, 
and evolve very fast, making their detection very difficult. 
The study of massive star formation also poses major theo- 
retical challenges because they begin burning their nuclear 
fuel and radiating prodigious amounts of energy while still 
accreting. If the formation of massive stars is similar to that 
of low mass stars (i.e. via accretion from the surrounding 
envelope), a mass accretion rate M of at least several or- 
ders of magnitude above the values appropriate for the low 
mass star formation is required (Tan & McKee 2002). At 
present, two theoretical scenarios are proposed to explain 
the formation of these stars: a monolithic collapse of turbu- 
lent gas on the scale of massive dense cores (Tan & McKee 
2002), which is a scaled- up version of the low- mass star for- 
mation picture, and a competitive one where accretion oc- 
curs inside the gravitational potential of a cluster-forming 
massive dense core (Bonnell & Bate 2006). This last model 
predicts that stars located near the center of the full grav- 
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itational potential accrete at much higher rates than do 
isolated stars. 

The detection of molecular outflows associated with 
both, high and low-mass young stellar objects, supports 
a scaled- up formation picture (Beuther 2002a). In this con- 
text, the study of massive molecular outflows and the asso- 
ciated driving source might contribute to discern which is 
the scenario that prevails in a given star forming region. 

It is well known that the generation of molecular out- 
flows during the formation of a high-mass star is a phe- 
nomenon that can take place even when the UCHII re- 
gion stage has been reached (Hunter et al. 1997; Qin et al. 
2008). In this work we report the study of the UCHII re- 
gion G045.47+0.05 (hereafter G45.47), through the analysis 
of its associated molecular gas and searching for evidence 
of molecular outflows and its associated driving source. 
G45.47 was first detected by Wood & Churchwell (1989) 
in radio continuum at 6 cm. This object is adjacent to the 
extensively studied UCHII region G45.45+0.06 (hereafter 
G45.45), which is part of a complex of five radio compact 
HII regions (Matthews et al. 1977; Giveon et al. 2005a; 
Giveon et al. 2005b). Such complex and G45.47 are em- 
bedded in the molecular cloud GRSMC G045.49+00.04 (at 
^lsr ~ 58 km s _1 ; Rathborne et al. 2009) and are located 
on the northern border of the more extended HII region 
named G45L in Paron et al. (2009). Based on HI absorp- 
tion profiles, Kuchar & Bania (1994) derived a kinematics 
distance of 8.3 kpc for the UCHII region G45.45. As G45.47 
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is part of the same complex, we adopt for this object the 
same distance. 

Caswell et al. (1995) detected a class II CH3OH maser 
emission at 6.6 GHz at Vlsr ^ 56 km s _1 towards G45.47. 
Given that the 6.6 GHz methanol maser is radiatively 
pumped by IR emission from the warm dust associated ex- 
clusively with massive young stellar object (MYSOs) 1 their 
detection is useful to study the kinematic of the gas and, in 
particular, to establish the systemic velocity (Cyganowski 
et al. 2009). 

Based on high resolution molecular line observations to- 
wards G45.47, Wilner et al. (1996) identified five HCO+ 
(1-0) clumps and suggested that G45.47 is in the early 
stages of the formation of an OB cluster. From an ammo- 
nia absorption study, Hofner et al. (1999) suggested the 
presence of a remnant molecular core infalling onto the 
UCHII region. Later, Cyganowski et al. (2008) identified an 
"Extended Green Object" (EGO) at the position of G45.47, 
the EGO G45.47+0.05. Their identification as EGOs comes 
from the common coding of the 4.5 jam band as green in the 
three-color composite Infrared Array Camera (IRAC; Fazio 
et al. 2004) images from the Spitzer Telescope. Extended 
4.5 fim emission is thought to evidence the presence of 
shocked molecular gas in protostellar outflows. The asso- 
ciation of EGOs with IRDCs and 6.7 GHz CH 3 OH maser 
suggests that EGOs trace the formation of massive proto- 
stars. According to Cyganowski et al. (2008), an EGOs is 
a MYSO with ongoing outflow activity and actively accret- 
ing. 

In summary, G45.47 is a rich and complex region to 
study the formation of a new generation of massive stars. 
In this paper we present a new study of the dense ambient 
medium where the UCHII region is evolving. We investigate 
the molecular gas through several molecular lines observed 
with the Atacama Submillimeter Telescope Experiment 
(ASTE; Chile) and characterize the central source based 
on infrared public data. 

2. Data 

The molecular line observations were carried out on June 
12 and 13, 2011 with the 10 m Atacama Submillimeter 
Telescope Experiment (ASTE; Ezawa et al. 2004). We used 
the CATS345 GHz band receiver, which is a two-single band 
SIS receiver remotely tunable in the LO frequency range of 
324-372 GHz. We simultaneously observed 12 CO J=3-2 
at 345.796 GHz and HCO+ J=4-3 at 356.734 GHz, map- 
ping a region of 2 / x2 / centered at RA = 19 h 14 m 25.6 s , 
dec.=+ll o 09'27.6"(J2000). We also observed 13 CO J=3-2 
at 330.588 GHz and CS J=7-6 at 342.883 GHz towards 
the same region. The mapping grid spacing was 20" in 
both cases and the integration time was 20 sec ( 12 CO 
and HCO + ) and 40 sec ( 13 CO and CS) per pointing. 
All the observations were performed in position switching 
mode. We verified that the off position (RA = 19 h 14 m 21.6 s , 
dec.=+10°59 / 4 // , J2000) was free of emission. We used the 
XF digital spectrometer with a bandwidth and spectral res- 
olution set to 128 MHz and 125 kHz, respectively. The 
velocity resolution was 0.11 km s _1 and the half-power 
beamwidth (HPBW) was 22", for all observed molecular 

1 We define MYSOs to be young stellar objects (YSOs) that 
will eventually become main-sequence O or early B type stars 
(M* > 8 M ). 



lines. The system temperature varied from T sys = 150 to 
200 K. The main beam efficiency was ?7 m b ~0.65. All the 
spectra were Hanning smoothed to improve the signal-to- 
noise ratio. The baseline fitting was carried out using sec- 
ond order polynomials for the 12 CO and 13 CO transitions 
and fourth order polynomials for the HCO + and CS tran- 
sitions. The polynomia was the same for all spectra of the 
map of a given transition. The resulting rms noise of the 
observations was about 0.2 K. The data were reduced with 
NEWSTAR 2 and the spectra processed using the XSpec 
software package 3 . 

The observations are complemented with near- and mid- 
IR data extracted from public databases and catalogues, 
which are described in the corresponding sections. 

3. Results and Discussion 

Figure l-(left) shows a composite Spitzer-IRAC three-color 
image (3.6 //m = blue, 4.5 //m = green, and 8 //m = 
red). The extended HII region G45L is centered at RA 
= 19 h 14 m 17 s , dec.=+ll o 07'48"(J2000) and is delimited by 
two arclike structures observed at 8 /im. The UCHII re- 
gion G45.47 (a.k.a. EGO G045.47+0.05) is the green source 
located at RA=19 h 14 m 25.6 s , dec.=+ll°09'27.6"(J2000). 
The white box indicates the region mapped with ASTE. 
A zoom up of this region (Fig. 1- (right)) shows that 
the UCHII region is embedded in the Spitzer dark cloud 
SDC G45.467+0.048 (Peretto & Fuller 2009). Infrared dark 
clouds (IRDCs) are dense molecular clouds which appear as 
extinction features against the bright mid-infrared Galactic 
background and have been suggested as the cold precursors 
to high- mass stars (Rathborne et al. 2006). 

3.1. The molecular gas 

In this section, we present the molecular results starting the 
description with the transitions that trace the most inner 
part of the clump related to G45.47, and moving to those 
that map the external layers which may give information 
about the dynamic effects occurring in the gas. 

Figure 2-a shows the CS J=7-6 spectra obtained to- 
wards the 2' x 2' region (white box in Fig. 1-left). The 
mapped area includes G45.47 and part of the nearby UCHII 
region G45.45, located at the (0, 0) and (-60, -20) off- 
set, respectively. As expected, the most prominent emission 
arises from these two position, since the detection of the CS 
J=7-6 transition reveals the presence of warm and dense 
gas. Figure 2-b shows an enlargement of the spectrum ob- 
served at (0, 0). The emission related to G45.47 shows a 
triple peak structure with velocity components centered at 
about 56 (the systemic velocity of the gas), 62, and 65 km 
s _1 (see Table 1) with a pronounced dip at about 59 and a 
less conspicuous one at ~64 km s _1 . The velocity compo- 
nent centered at 56 km s _1 is weak, with a peak about 3<r of 
the rms noise level. Since the probability of superposition 
along the line of sight of more than one component in the 
CS J =7-6 line is very low, we suggest that the dips reveal 

2 Reduction software based on AIPS developed at NRAO, ex- 
tended to treat single dish data with a graphical user interface 
(GUI). 

3 XSpec is a spectral line reduction package for astronomy 
which has been developed by Per Bergman at Onsala Space 
Observatory 
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Fig. 1. Spitzer-IRAC three color images (3.6 /xm = blue, 4.5 /im = green and 8 fim. — red). The white box in the left panel indicates 
the region mapped with ASTE. A close up view of this region is shown in the right panel. Red, green and blue scales go from 30 
to 500, from 2 to 110, and from 1 to 80 MJy/sr, respectively. 



self-absorption effects in the gas, consistent with an opti- 
cally thick transition. Self-absorption features demonstrate 
the existence of a density gradient in the clump (Hiramatsu 
et al. 2007) and is a common feature of optically thick lines 
in the direction of embedded young stars. The depression 
discloses the presence of relatively cold, foreground gas that 
absorbs photons from warmer material behind it. The CS 
J=7-6 profile is asymmetric with respect to the dip near 
v=59 km s _1 as the redshifted component is clearly brighter 
than the blueshifted one, suggesting that the molecular gas 
is expanding. This is because in an expanding cloud a line 
emission is composed by red and blueshifted photons, the 
redshifted photons will encounter fewer absorbing material 
(which is expanding outward) than would blueshifted pho- 
tons and hence have greater probabilities of escaping (e.g. 
Leung 1978; Zhou 1992; Lehtinen 1997). 

Figure 3 shows the velocity channel maps of the CS 
J=7-6 emission averaged every 1 km s _1 . It can be noticed 
the presence of molecular gas associated with G45.47 (red 
cross) in the velocity range going from ~ 55 to ~65 km 
s _1 . The other most conspicuous molecular condensation, 
partially observed, corresponds to the gas related to G45.45 
and has a central velocity of ~ 59 km s _1 . 

Figure 4- a shows the HCO + J=4-3 spectra obtained 
towards the same region. As in the case of the CS J=7-6 
spectrum it can be seen that the HCO + emission towards 
(0, 0) can be deconstructed in three velocity components, 
a very weak one centered at about 56 km s _1 (intensity of 
about 2<r), plus two brighter ones at 61 and 65 km s _1 , 
with only one dip centered at ~ 63 km s -1 . 

From Fig. 4- a it can be noticed that the HCO + spec- 
trum on position (—60, —20), that is at the position of 
the UCHII region G45.45, has also a dip at about 60 km 
s _1 similarly to what it is observed in the CS J=7-6 line 
towards the same position. 

The HCO + spectra also show evidence of emission 
near position (60, —40). This molecular feature is in po- 



sitional coincidence with an infrared (IR) source located 
at RA=19 h 14 m 27.7 s , dec.=+ll°08 / 33 // (J2000) (see Fig 1). 
As this HCO + spectrum has the same kinematic structure 
than the spectrum towards (—60, —20) including the pres- 
ence of a dip at the same velocity of about 60 km s -1 , we 
can conclude that this IR source is embedded in a molecular 
filament that must be connected with G45.45. 

Figure 5 shows the velocity channel maps of the HCO + 
J=4-3 emission averaged every 1 km s -1 . The HCO + J=4- 
3 emission related to G45.47 is visible between ~ 55 and ~ 
68 km s -1 , while the emission associated with G45.45 goes 
from ~ 54 to ~ 63 km s _1 , respectively. Between ~ 56 and 
~ 60 km s _1 part of the molecular condensation related to 
the IR source mentioned above can be appreciated. 

Figure 6 shows the 13 CO J=3-2 spectra towards the 
same region. The 13 CO emission towards (0, 0) can be de- 
composed in three velocity components centered at about 
56 (the systemic velocity of the gas), 60, and 63 km s -1 , 
while only one dip centered at ~ 58 km s _1 can be appre- 
ciated. The depression at 58 km s _1 can be observed in 
almost all the 13 CO J=3-2 spectra in the region. Figure 7 
shows the velocity channel maps of the 13 CO J=3— 2 
emission averaged every 1 km s _1 . 

The 12 CO J=3-2 spectra (Fig. 8) exhibit a more compli- 
cated behaviour. The profile towards the position (0, 0) has 
a triple peak structure with components centered at about 
55, 60, and 64 km s _1 and two dips at about 58 and 62 km 
s -1 . This spectrum is broadened, suggesting the presence 
of outflowing activity in the region with the blue wing cen- 
tered near the position (20, 20) and the red wing around 
the (0, -40) offset. 

In Figure 9 we show the velocity channel maps of the 
12 CO J=3-2 emission averaged every 2.1 km s _1 . Among 
all the observed molecular condensations, we draw the at- 
tention onto the clumps related to G45.47 and G45.45. 
The molecular clump related to G45.47 is observed be- 
tween ~ 62 and 66 km s _1 and is seen slightly shifted to 
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Fig. 2. a)CS J=7-6 spectra obtained towards the 2' x 2' region (white box in Fig. 1-left) mapped with ASTE. b) Hanning 
smoothed profile of the CS J=7-6 line towards G45.47 at the position (0, 0). The dashed line indicates the systemic velocity of 
the gas. The spectrum was deconstructed using three Gaussians, which are shown in red. 




Fig. 3. Velocity channel maps of the CS J=7-6 emission averaged every 1 km s 1 . Greyscale goes from 0.15 (about 3a rms noise 
of an averaged channel map) up to 1.8 K. Contours levels are at 0.15, 0.3, 0.45 and 0.6 K. The red cross represents the position 
of G45.47. 



the northwest between ~ 58 and 62 km s 1 . The clump 
associated with G45.45 (partially observed) is centered at 
RA=19 h 14 m 22 s , dec.=+ll o 09'20"(J2000) in the velocity 
interval going from 47 to 63 km s _1 . Although both molecu- 
lar condensations have associated different velocity ranges, 
they are connected through the extended emission shown in 
the velocity interval going from 50 to 65 km s _1 . The molec- 
ular gas related to the spectral wings appears as two con- 
spicuous molecular features centered at RA=19 h 14 m 27 s , 
dec.=+ll°09'45"(J2000) between 35 and 53 km s" 1 and 
at RA=19 h 14 m 26 s , dec.=+ll o 08'45"(J2000) between 65 
and 76 km s _1 . These features will be further discussed 
in Section 3.3. 



The four transitions have, within errors, the same main 
velocity components at about 55, 60, and 64 km s _1 to- 
wards the position (0, 0). By the other hand, the velocity 
components observed at about 51 and 68 km s _1 in the 
12 CO J=3-2 spectrum are not detected in the other three 
lines. Table 1 lists the emission peaks parameters derived 
from a Gaussian fitting for the four molecular transitions 
on the position (0, 0). T m 5 represents the peak brightness 
temperature and Vlsr the central velocity referred to the 
Local Standard of Rest. Errors are formal la value for the 
model of Gaussian line shape. All the spectra towards the 
(0,0) position have the same self- absorption dip at about 
58-59 km s _1 , which correspond to the central velocity of 
the molecular cloud GRSMC 045.49+00.04 where G45.47 is 
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Fig. 4. a) HCO + J=4-3 spectra obtained towards the 2' x 2' region (white box in Fig. 1-left) mapped with ASTE. b) Hanning 
smoothed profile of the HCO + J=4-3 line towards the position (0, 0) where it is G45.47. The dashed line indicates the systemic 
velocity of the gas. The spectrum was deconstructed using three Gaussians, which are shown in red. 




Fig. 5. Velocity channel maps of the HCO + J= 4-3 emission averaged every 1 km s -1 . Greyscale goes from 0.15 (about 3a rms 
noise of an averaged channel map) to 2.5 K. Contours levels are at 0.15, 0.5, 0.8, 1.3, and 1.9 K. The red cross represents the 
position of G45.47. 



embedded. Even more, this spectral feature is observed in 
all transitions towards the region G45.45 and the IR source 
located at RA=19 h 14 m 27.7 s , dec.=+ll°08'33"(J2000) (see 
Fig 1). 

3.2. Column density and mass estimates of the molecular 
clump associated with G45.47 

We estimate the 13 CO J=3-2 opacity, T13, based on the 
following equation: 

7 / T peafc ( 13 CO) ^ 



where we consider T pea k from the position (0, 0). We ob- 
tain, T13 ~ 1.9 which reveals that the 13 CO J=3-2 emission 
is optically thick towards G45.47, in agreement with the ob- 
served profile towards (0,0) offset (see Fig. 6). 

The excitation temperature, T eXl of the 13 CO J=3-2 
line was estimated from: 



Tpeak k \e h »/ kT ™ - 1 e h »/ kT BG - 1 J X ^ 6 Tl3 ) 

(2) 

where for this line hv/k = 15.87. Assuming Tbg — 2.7 K, 
and considering the peak brightness temperature for the 
13 CO J=3-2 at (0, 0) offset, 13 T peak = 11.15 K, we derive 
a T ex - 20 K for the 13 CO J=3-2 line. 
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Fig. 6. 13 CO J=3-2 spectra obtained towards the 2' x 2' region (white box in Fig. 1-left) mapped with ASTE. b) Hanning 
smoothed profile of the 13 CO J=3-2 line towards G45.47. The dashed line indicates the systemic velocity of the gas. The spectrum 
was deconstructed using three Gaussians, which are shown in red. 




Fig. 7. Velocity channel maps of the 13 CO J=3-2 emission averaged every 1 km s 1 . Greyscale goes from 0.15 up to 9 K. Contours 
levels are at 0.6, 2, 4, 5, 6, 7, 8, and 9 K. The red cross represents the position of G45.47. 



Finally, using the RADEX 4 code (van der Tak et al. 
2007) we derive the 13 CO J=3-2 column density and the 
H2 volume density. The RADEX model uses the mean es- 
cape probability approximation for the radiative transfer 
equation. 

Adopting T ex ~ 20 K, r 13 ~ 1.9, and 13 T peak = 
11.15 K, we obtain a 13 CO column density N( 13 CO) ~ 
2.8 x 10 17 cm -2 and n(H2) ~ 10 5 cm -3 . Considering an 
abundance ratio of [H 2 ]/[ 13 CO] = 77xl0 4 (Wilson & Rood 
1994) we estimate the H2 column density, N(H2) ~ 2.1 x 
10 23 cm -2 . Finally, using the relation M = /ira#<i 2 ^N(H2), 
where fi is the mean molecular weight per H2 molecule 

4 http:/ /www. sron.rug.nl/~vdtak/radex/radex.php 



(fi ~ 2.72), rriH the hydrogen atomic mass, d the distance, 
and ft the solid angle subtended by the structure, we cal- 
culate the total mass of the clump in M ~ 10 4 M©. 

As an independent estimate, we derive the beam- 
averaged gas column density, the mass, and the volume 
density of the clump based on the dust continuum emission. 
In particular, we use the integrated flux of the continuum 
emission at 1.1 mm as obtained from The Bolocam Galactic 
Plane Survey II Catalog (BGPS II; Rosolowsky et al. 2010). 
The 1.1 mm continuum emission is mostly originated in 
optically thin dust (Hildebrand 1983). Following Beuther 
et al. (2002a) and Hildebrand (1983) we calculate the mass 
and the gas column density of the clump using: 
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Fig. 8. a) 12 CO 3=3-2 spectra obtained towards the 2' x 2' region (white box in Fig. 1-left) mapped with ASTE. b) Hanning 
smoothed profile of the 12 CO J=3-2 line towards G45.47. The dashed line indicates the systemic velocity of the gas. The spectrum 
was deconstructed using five Gaussians, which are shown in red. 
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Fig. 9. Velocity channel maps of the 12 CO J=3-2 emission averaged every 2.1 km s _1 . The given velocities correspond to the 
higher velocity of each interval. Greyscale goes up to 23 K. Contours are above the 5cr of the rms noise level. The cross indicates 
the position of G45.47. 
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Fig. 10. Comparison between the 12 CO J=3-2 and the HCO + J=4-3 spectra on the positions (20, 20) and (0, -40). The dashed 
line indicates the systemic velocity. The velocity intervals corresponding to the blue and red wings are shown. The four spectra 
were Hanning smoothed. 



Table 1. Emission peaks parameters derived from a Gaussian 
fitting for the four molecular transitions spectra on the position 
(0, 0). 



Transition 


Vlsr [km s -1 ] 


T mb [K] 


Av [km s -1 ] 


CS J=7-6 


56.3±0.9 


0.6±0.2 


3.1±0.6 




62.0±0.8 


1.8±0.1 


3.1±0.7 




64.8±1.1 


1.3±0.1 


2.5±0.9 


HCO+ J=4-3 


55.8±1.2 


0.4±0.2 


5.9±0.5 




60.8±1.4 


4.4±0.8 


3.8±0.6 




64.7±1.7 


5.1±0.6 


3.3±0.7 


L6 CO J=3-2 


55.9±0.6 


5.2±0.4 


2.9±0.4 




60.5±0.5 


11.1±0.8 


3.3±0.6 




63.3±0.4 


6.4±0.7 


2.6±0.7 


A2 C0 J=3-2 


50.7±0.3 


3.0±0.2 


9.0±0.3 




55.2±0.1 


13.1±1.1 


3.3±0.4 




59.8±0.2 


10.6±0.8 


3.6±0.6 




64.4±0.6 


10.6±0.7 


3.3±0.5 




68.1±0.4 


1.3±0.2 


4.4±0.4 



_ 1.3 x 10" 3 a p R F v 
9as ~ JATdust) 0.1/mi3gcm- 3 100 

and 



7.8 x 1Q 1Q a _p R_F^_ 

gaS ~ Ju(T dust )n b 0.1/im 3gcm- 3 100 

where J v (T dust ) = [exp(his / kT dust ) - l] -1 and Q b ,a,p,R, 
and (3 are the beam solid angle, grain size, grain mass 



density, gas-to-dust ratio, and grain emissivity index for 
which we used the values of (33") 2 in radians, 0.1 /im, 
3 g cm -3 , 100, and 2, respectively (Hunter 1997, Hunter 
et al. 2000, and Molinari et al. 2000). Based on the work 
of Sridharan et al. (2002) who derived dust temperatures 
ranging between 30 and 60 K for a sample of several mas- 
sive star forming regions, we adopt T dus t = 45 K. For a 
distance d = 8.3 kpc and an integrated flux intensity F v 
= 5.18 Jy at 1 mm (Rosolowsky et al. 2010) we obtain 
Ngas - 4 x 10 23 cm -2 , M gas - 8520 M , and a volume 
density, n(H.2) ~ 1.4 x 10 5 cm -3 . These values are in good 
agreement with those derived from the 13 CO J=3-2 line 
using RADEX. 

3.3. Molecular outflows associated with G45.47 

As discussed in Section 3.1, the presence of spectral wings 
in the 12 CO J=3-2 spectrum obtained towards G45.47 is 
a strong indicator that molecular outflow activity is taking 
place in the region. To characterize the associated molecular 
outflows it is first necessary to separate the outflowing gas 
from the molecular material of the clump, identifying the 
velocity ranges related to each structure. 

We consider two independent methods to determine 
these velocity ranges. The first method consists in the iden- 
tification of the blue and red spectral wings based on a com- 
parison between the 12 CO J=3-2 and HCO + J=4-3 spec- 
tra. Figure 10 shows the comparison between both spectra 
on positions (20, 20) and (0, —40) where the blue and red 
wings are largest. Considering emission up to about 2 a of 
rms noise level, it can be noticed a blue and a red wing in 
the 12 CO spectra between ~ 34 and ~ 54 and between ~ 
64 and ~ 75 km s _1 , respectively. 

The second method to identify the molecular outflows is 
based on the inspection of the 12 CO J=3-2 data cube, chan- 
nel by channel, trying to spatially separate both outflows's 
lobes. In Section 3.1 we mentioned two molecular features 
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in the 12 CO J=3-2 emission centered at RA=19 h 14 m 27 s , 
dec.=+ll o 09'45"(J2000) between - 34 and - 54 km 
s- 1 and at RA=19 h 14 m 26 s , dec.=+ll o 08'45"(J2000) be- 
tween ~ 64 and ~ 75 km s _1 which we identify as the out- 
flow's lobes related to G45.47 (see Fig. 9). These molecular 
features can be identified as the spectral wings detected in 
the 12 CO J=3-2 spectra through the first method. Figure 
11 shows a Spitzer-IRAC three color image (3.6 jiva = blue, 
4.5 /im = green and 8 fim = red) of G45.47. The blue and 
red contours represent the 12 CO J=3-2 emission integrated 
from ~ 34 to 54 km s _1 (blue lobe) and from ~ 64 to 75 km 
s _1 (red lobe), respectively. 




Fig. 11. Spitzer-IRAC three color image (3.6 /im = blue, 4.5 /mi 
= green and 8 /am = red) of G45.47. The blue and red contours 
represent the 12 CO J=3-2 emission integrated from 34 to 54 km 



* (red lobe), respec- 
600, and 670 K km 



s (blue lobe), and from 64 to 75 km s 
tively. The blue contours are at 400, 500, 
s" 1 and the red ones are at 140, 160, 180, and 200 K km s -1 . 

Following Scoville et al. (1986) we estimate the 12 CO 
optical depth, 712, of the gas in the molecular outflows us- 
ing: 



exp{-ri2) 



■ exp{—T\<ijX) 



(5) 



where X=[ 12 CO]/[ 13 CO] is the isotope abundance ratio. 
Using the relation X = 6.21 x D GC + 18.71 (Milam et al. 
2005), where Dgc — 6.1 kpc is the distance between G45.47 
and the Galactic Center, we obtain an isotope abundance 
ratio X— 57 for this region. We adopt a constant X 
throughout the outflows (Cabrit & Bertout 1990). Since 
we could not identify any blue wing in the 13 CO spectrum 
on position (20, 20), the 12 T m 5/ 13 T m 5 ratio for both wings 
was estimated considering only the red one. We obtain a 
12 T m 6/ 13 T m 6 ratio of about 4.2 and derive a 12 CO optical 
depth of r\2 ~ 15, which will be adopted for both spectral 
wings in further calculations. 

We can now estimate the 12 CO column density of both 
outflow lobes from (see e.g. Buckle et al. 2010): 



N( 12 CO) = 7.96 x 10 13 e^ 



13 ^e T ex + 0.92 
l-exp(=^) 



J ri 2 dv (6) 



Taking into account that the 12 CO J=3-2 is an optically 
thick transition (r > 1), the integral can be approximated 
by: 



/ 



Tl2dv 



Tl2 



J(T e; 



J(T BG ) 1 - 



l Tm h dv (7) 



By integrating the 12 CO emission from 34 to 54 km s 1 and 
from 64 to 75kms _1 we obtain N (CO) Uue ~ 5.2xl0 17 cm~ 2 
and N(CO) re( i ~ 2 x 10 1T cm -2 , respectively. The mass 
of each wing can be derived from the relation M = 
/im // d 2 ^X(CO)- 1 N(CO), where X(CO) is the 12 CO rel- 
ative abundance to H 2 (X(CO)~ 7.4 x 10" 5 ). We obtain 
M biue ~ 300 M and M red ~ 120 M , yielding a total 
outflow mass of about 420 M . In this way, the estimated 
total outflows mass represents about 4% of the clump mass, 
which is in agreement with the results of Beuther (2002b) 
who established that approximately 4% of the clump gas is 
entrained in the molecular outflows. 

The momentum and the kinetic energy of the wings can 
be derived from P = MV C and E k = O.SMV^ 2 , where V c is a 
characteristic velocity estimated as the difference between 
the maximum velocity of detectable 12 CO emission in the 
wing and the systemic velocity of the gas (~ 56 km s _1 ). 
Taking into account a V blue ~ 22 km s _1 and a V^ ed ~ 
19 km s _1 , we obtain Puue = 6.6 x 10 3 M km s _1 , P re d = 
2.3 x 10 3 M km s"\ E b k lue = 1.4 x 10 48 erg and E r k ed = 
4.2 x 10 4T erg. Comparing our results with the work of Wu 
et al. (2004) who based on compiled data concluded that the 
mass, momentum, and energy of molecular outflows range 
from 10~ 3 to 10 3 M , 10~ 3 to 10 4 M km s -1 and 10 38 to 
10 48 erg, respectively, we can infer that in this work we are 
dealing with massive and energetic outflows. 

Following Curtis et al. (2010) we define the dynamical 
time for the blue and red wings, tdyn-> as the time for the 
bow shock travelling at the maximum velocity in the flow, 
V c , to travel the projected lobe length, Riobe' 



tdy 



R, 



lobe 



Vr 



(8) 



with Rf H = 1.8 pc and R t 



![^ e = 1.6 pc. These values were 
obtained by inspecting the blue and red lobe sizes in Fig. 11 
and considering a distance of 8.3 kpc. We obtain a similar 
dynamical time of both lobes, td yn — 0.8 x 10 5 yr. According 
to Beuther et al. (2002b), flow ages are good estimates of 
protostar lifetimes. We can then assume that the protostar 
responsible of the outflow activity has an age of the order of 
10 5 yr. This is consistent with the detection of ionized gas 
in the region, since to be able to ionize its surroundings the 
age of the protostar should be at least 10 5 yr (Sridharan 
et al. 2002). 



3.4. HCO+ as tracer of molecular outflows 

As an independent test of out findings, we performed a 
study of the HCO + abundance towards G45.47. HCO + is 
believed to be the dominant ionized species in dense dark 
clouds (Rawlings et al. 2000; Dalgarno & Lepp 1984), being 
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very important for the ion-neutral chemistry. The ioniza- 
tion fraction of molecular clouds is a relevant parameter 
to study the cloud chemistry and dynamics, and hence the 
star forming processes. Star formation occurs in the inte- 
rior of dense cores, regions of high extinction where self- 
shielding prevents the UV photoionization of H2. Thus, 
cosmic ray ionization is believed to dominate photoioniza- 
tion in dense cores (McKee 1989). However, when YSOs 
are formed within the cores, the shocks and outflows pro- 
duce deep changes in the chemistry and photoionization 
processes are very likely. Indeed, in star forming regions 
occur a molecular enrichment due to the desorption of 
molecular-rich ice mantles, followed by photochemical pro- 
cessing by shock-generated radiation fields. In particular, 
towards YSOs which are driving outflows, it is expectable 
an enhancement in the HCO + abundance (Rawlings et al. 
2000, 2004). In what follows, we estimate, using a simple 
chemical network, the HCO + abundance that would be pro- 
duced by an standard cosmic ray ionization rate in the G45 
clump in order to compare with the abundance obtained 
from our observations. 

The starting chemical reaction to form HCO + is the 
production of H3" , mainly formed in the interaction between 
the cosmic rays (c.r.) and the molecular gas (e.g. Oka 2006): 



H 2 + c.r. -> + e 
H+ + H 2 -> H+ + H. 

The rate equation of this reaction is 



Ch 2 n(R 2 ) 



dt 



(9) 



(10) 



where (h 2 is the cosmic ray ionization rate and n the density 
of the molecular species. By considering that the HCO + is 
mainly formed by the reaction of H3" with CO: 



H+ + CO -> HCO+ + H 2 



(11) 



and destroyed through recombination with electrons, 

HCO + +e- -^CO + H, (12) 

the rate equations can be equated, leading to 

^hco+ n(H+) n(CO) = k CO n(HCO+) n(e), (13) 

where /c H co+ an d kco are the coefficient rates and n(e) 
the electron density. Assuming that the main destruction 
mechanism for Hjj" is through the formation of HCO + , the 
rate of destruction of H3" is therefore equal to the formation 
rate of HCO + ; this implies that the left side of equation (10) 
is equal to right side of equation (13). Then, it is possible 
to write an expression for the cosmic ionization rate as a 
function of the molecular densities: 



Ch 2 = &co 



n(HCO+) n(e) 



(14) 



The rate coefficient fcco? extracted from the UMIST 
database (Woodall et al. 2007), is: k C o = 2.4 x 
10- 7 (T/300K)-° 69 . Equation (14) can be approximated 
using the column densities, leading: 



Ch 2 — ^co 



N(HCO+) X(e) n(H 2 ) 
N(H 2 ) ' 



(15) 



where X(e) is the electron abundance. Using this equa- 
tion and assuming typical values for the cosmic ionization 
rate and electron abundance of £h 2 = (1 — 5) x 10 -17 s _1 
(Dalgarno 2006) and X(e) ~ 10~ 7 (Bergin et al. 1999), 
respectively, we derive the HCO + column density. We use 
N(H 2 ) = 2 x 10 23 cm" 2 and n(H 2 ) = 1 x 10 5 cm" 3 , as de- 
rived above, and we assume T = 20 K. Finally, we obtain 
that the N(HCO+) should be in the range (1.5-6.5) x 10 14 
cm -2 , leading to an abundance X(HCO + ) in between 
6.5 x 10- 10 and 3.2 x lO" 9 . 

On the other hand, we analyze the central HCO + J=4- 
3 spectrum and use the RADEX code (van der Tak et al. 
2007) to derive its column density. Assuming background 
and kinetic temperatures of 2.73 K and 20 K, respec- 
tively, and the same H 2 density used above, we varied 
the HCO + column density until obtaining a good fit for 
the observed peak temperature. The best fit was obtained 
with N(HCO + )~ 1.5 x 10 15 cm -2 , leading an abundance 
of X(HCO+) - 7.5 x 10~ 9 . This value, more than twice 
greater than the lager value obtained above, might indicate 
that the cosmic ray ionization is insufficient to produce the 
observed HCO + abundance through the described chemical 
network. Therefore, following Rawlings et al. (2000, 2004), 
we suggest that the observed HCO + abundance must be 
mostly produced by the outflowing activity in the region. 
In spite of the uncertainties in the calculations, this is an 
independent proof pointing to support an scenario with out- 
flows in the region. 

3.5. Looking for the driving source of the molecular outflows 

Wilner et al. (1996) reported 5" resolution observations of 
the HCO + J=l-0 transition towards G45.47 obtained with 
the OVRO millimeter array. They identified at least five 
HCO + J=l-0 clumps and suggested that the formation of 
an OB cluster would be taking place in the region. However 
the authors did not carry out any search of infrared point 
sources in the region to confirm that hypothesis. In this 
context, we wonder if the massive molecular outflows ob- 
served towards G45.47 were originated by one or more stars. 
Based on a near-infrared analysis, we searched for the driv- 
ing source candidates of the massive molecular outflows. 

Figure 12 shows the radio continuum emission at 6 cm 
obtained from The Multi-Array Galactic Plane Imaging 
Survey (MAGPIS; White et al. 2005) of the area mapped 
using ASTE. The green contours represent the HCO + J=4- 
3 emission distribution integrated between 58 and 67 km 
s _1 . The black contours were taken from the paper of 
Wilner et al. (1996) and represent the HCO + J=l-0 emis- 
sion distribution integrated in the same velocity interval. 
The region observed by Wilner et al. (1996) is indicated by 
the dashed box. The five clumps reported by the authors 
have been labeled. The positional coincidence among the 
center of the HCO + J=4-3 clump, the clump 3 of HCO + 
J=l-0 and a radio continuum source, which has been identi- 
fied as the UCHII region G45.47, is striking. The red circles 
indicate the location of the YSO candidates in the observed 
region (see red circles in Fig. 13). Among the five YSO can- 
didates projected onto the HCO + J=4-3 emission, 2MASS 
19142564+1109283 is the only one having a positional co- 
incidence with a clump of HCO + J=l-0 (clump 3). We do 
not find any embedded infrared sources in the others four 
clumps, suggesting that these clumps could be in a prestel- 
lar core stage or they might be tracing the origin of the 
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molecular outflows. We suggest that the most likely candi- 
date to be the driving source of the molecular outflows is 
2MASS 19142564+1109283. 
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Fig. 12. Radio continuum emission at 6 cm (MAGPIS; White 
et al. 2005) of the area mapped using ASTE. The green contours 
represent the HCO + J=4-3 emission distribution integrated be- 
tween 58 and 67 km s _1 . The black contours were taken from 
the paper of Wilner et al. (1996) and represent the HCO + J=l- 
emission distribution integrated from 58 and 67 km s _1 . The 
box indicates the area observed by Wilner et al. (1996). The 
numbers show the five HCO + J=l-0 clumps. The red circles 
represent the YSO candidates shown in Fig. 13. 
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Fig. 13. Color-color diagram of the 2MASS infrared sources 
in the vicinity of G45.47. The two solid curves represent the 
location of the main sequence (thin line) and the giant stars 
(thick line) derived from Bessell & Brett (1988). The parallel 
dashed lines are reddening vectors with the crosses placed at in- 
tervals corresponding to five magnitudes of visual extinction. We 
have assumed the interstellar reddening law of Rieke & Lebofsky 
(1985) (Aj/Av =0.282; A H /A V =0.175 and A K /A V =0.112). 
The sources reddened by circumstellar dust are indicated as red 
circles. 



ified as upper limit do not contribute to N). Hereafter, we 
refer to models satisfying Eq. 16 as "selected models". 



Figure 13 shows a color-color diagram (CCD) including 
all the 2MASS point sources located in the observed region. 
The red circles represent the sources with infrared excess 
(YSO candidates) which have been shown in Figure 12. The 
reddenest object is 2MASS 19142564+1109283. The blue 
circles represent main sequence or giant star candidates. 

To characterize the infrared source 2MASS 
19142564+1109283 we perform a fitting of its spec- 
tral energy distribution (SED) using the tool developed 
by Robitaille et al. (2007) 5 . We adopt an interstellar 
extinction in the line of sight, A v , between 5 and 17 
magnitudes. The range of A v was chosen by inspecting the 
location in a [H-K] vs [J-H] CCD of the 2MASS sources in 
the region (see Fig. 13). To construct the SED we consider 
the fluxes at the JHK 2MASS bands, Spitzer-IRAC 5.8 
and 8 fim bands, WISE 3.4, 4.6, 12 and 22 /im bands, 
MSX 14 and 21 /mi bands, and SCUBA 850 /im band. 
The fluxes of the datasets having lower angular resolution 
(MSX and SCUBA) were considered as upper limits. 
Figure 14 shows the best fitting SEDs models for 2MASS 
19142564+1109283. We select models that satisfies the 
condition: 

X 2 - xlst < ™, (16) 

where xtest 1S tne minimum value of the x 2 among all mod- 
els, and TV is the number of input data fluxes (fluxes spec- 

5 htt p : / /caravan . astro . wise . edu/ prot ost ars / 
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Fig. 14. Best fitting SED models for G45.47. The fill circles 
indicate the measured fluxes at the JHK 2MASS bands, Spitzer- 
IRAC bands at 5.8 and 8.0 /mi, and the WISE bands at 3.4, 
4.6, 12 and 22 /im. Triangles indicate fluxes considered as upper 
limits, as MSX bands at 14 and 21 /im and the SCUBA band at 
850 /im. Black and gray solid curves represent the best-fit model 
and the subsequent good fittings, respectively. The dashed line 
shows the best-fit model for a central source contribution in 
absence of circumstellar dust. 

From the SED fitting several model parameters can be 
obtained and, as expected, some can be better constrained 
than others. Taking into account the values yielding by the 
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"selected models" for several parameters, we derive the fol- 
lowing main results: 

— the total luminosity distribution of the source has an 
average value of 2 x 10 4 L with a spread between 8 x 10 3 
and 5 x 10 4 L©. 

— embedded in the molecular clump there is a massive 
protostar of about 15 M (early B-type star) accreting 
material from its envelope at large rates, M env ^ 2 x 
l(r 4 M yr- 1 . 

— despite the large spread in the central source's age dis- 
tribution that goes from 2 x 10 3 to 10 6 yr, it peaks at 
about 3 x 10 5 yr, as expected for protostars that has be- 
gun to ionize their surroundings. Besides, this result is 
in agreement with the dynamical time, td yn = 1 x 10 5 yr, 
derived in Section 3.3. 



4. Summary 

We carried out a study of the UCHII region G045.47+0.05 
and its surroundings based on molecular lines observations 
and public infrared data. We find a molecular gas conden- 
sation associated with G045.47+0.05. The detection of the 
CS J=7-6 transition reveals that the UCHII region is still 
embedded in warm and dense molecular gas. The CS spec- 
trum obtained towards the position (0, 0) has a pronounced 
self- absorption dip at about 59 km s _1 with the red velocity 
component more intense than the blue one. This asymme- 
try in the CS profile suggests an expansion of the molecular 
gas. 

All molecular transitions have the same main three ve- 
locity components with a self-absorption dip at about 58- 
59 km s _1 which correspond to the central velocity of the 
molecular cloud GRSMC G045.49+00.04 where G45.47 is 
embedded. 

Based on the ratio between the 12 CO and 13 CO J=3-2 
transitions we estimate the 13 CO opacity, r±s, towards the 
center of the clump in about 1.9, revealing that this line is 
optically thick. This result is in agreement with the struc- 
ture exhibited by the 13 CO spectrum towards the center of 
G045.47+0.05. 

Using the RADEX code we derive a 13 CO J=3-2 column 
density and H2 volume density of about 2.8 x 10 17 cm -2 
and 10 5 cm -3 , respectively. The H2 column density and the 
total mass of the clump were estimated in 2.1 x 10 23 cm -2 
and M ~ 10 4 M , respectively. 

From an independent estimate based on the dust con- 
tinuum emission at 1 mm we derive N gas ~ 4 x 10 23 cm -2 , 
M gas ~ 8520 M , and a volume density, n(H 2 ) ~ 1.4 x 
10 5 cm -3 , in good agreement with the values derived from 
the molecular lines. 

From the analysis of the 12 CO J=3-2 and HCO+ tran- 
sitions we report the presence of molecular outflows related 
to the UCHII region G045.47+0.05. We estimate the blue 
and red lobes masses in about 300 and 120 M , respectively. 
The total mass of the outflows represents about 4% of the 
clump mass. The dynamical time of both lobes, good indica- 
tor of the protostar age, was estimated in about 8 x 10 4 yr. 
This agrees with the presence of ionized gas around the 
central object. 

Based on infrared data we searched for the protostar (s) 
that must be driving the outflows. We find that the source 
2MASS 19142564+1109283 is a YSO candidate located 



onto the center of the clump in coincidence with the ra- 
dio continuum emission associated with the UCHII region. 
Besides this infrared source is the only YSO candidate re- 
lated to one of the five HCO + clump detected by Wilner 
et al. (1996). We did not detected any YSO candidate as- 
sociated with the other four HCO + clumps. This results 
suggests that the massive molecular outflows are generated 
by the protostar 2MASS 19142564+1109283. 

Finally, from a spectral energy distribution analysis of 
the source we derive a total luminosity, a protostar's mass, 
and protostar's age of about 2 x 1O 4 L , 15 M (early B-type 
star), and 3 x 10 5 yr, respectively. 

Massive molecular outflows give, on large scales, a good 
piece of information about the physical processes taking 
place in the innermost parts of the star- forming cores. Their 
detection give indirect evidence that massive star formation 
is a scaled-up version of low mass star formation. The iden- 
tification of a massive protostar that would be generating 
the massive outflows reinforces this scenario. 
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